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Matched filtering is a popular data analysis framework used to search for gravitational wave 
signals emitted by compact object binaries. The templates used in matched filtering searches are 
constructed predominantly from the quadrupolar mode because this mode is the energetically most 
dominant channel. However, for highly precessing binaries or binaries with moderately large mass 
ratios, significant power is also carried by higher-order modes. We investigate how the inclusion of 
higher modes in the templates increases the prospects for detecting gravitational waves. Specifically, 
we use numerical relativity waveforms from the late inspiral and coalescence of binary black holes 
to identify mode hierarchies that cover the sky of binary orientations. We show that the ordering in 
these hierarchies depends on the characteristics of the binary system and the mode strengths. Our 
study demonstrates that detecting moderately high precessing or unequal mass binaries requires the 
inclusion of higher modes in the templates banks. 

PACS numbers: 



I. INTRODUCTION 

Binary systems with black hole (BH) components will 
provide us with one of the strongest sources of gravita- 
tional waves (GWs) soon to be detected by ground-based 
interferometers pQ. GW observations are accompanied 
by dazzling engineering, theoretical, and data analysis 
challenges. One challenge that requires the close cooper- 
ation between source model and data analysis researchers 
is digging out the GW signals from the noisy data. For 
compact object binaries involving neutron stars and BHs, 
one of the most popular frameworks to analyze the data 
is matched filtering, which in the broad sense consists 
of correlating a known signal or template with the data 
in hopes of identifying the presence of the template in 
the data. Not surprisingly, the success of matched filter- 
ing depends mainly on the "quality" of the template, in 
other words, the degree to which the template captures 
the correct physics and matches the characteristics of the 
source. Even with the belief that we possess the correct 
theory of gravity, Einstein's theory of general relativity, 
constructing templates is a laborious and arduous task. 
A family of binary black hole (BBH) templates needs 
to sample a 15-dimcnsional parameter space (BH masses 
and spins, binary eccentricity, orientation vector, sky po- 
sition and distance). Furthermore, if the binary merges 
in the sweet spot of the detector, building templates re- 
quires numerical relativity (NR) simulations that capture 
dynamically strong gravity, thus increasing the compu- 
tational cost per template significantly. 

Without accounting for effects from the orientation of 
the source, the gravitational radiation from a binary sys- 
tem is dominated by its quadrupolar (2,2) mode compo- 
nent. It should not be surprising then that current data 
analysis pipelines are predominantly using this mode as 
a template. An example of the dominance of the (2,2) 
mode is depicted in Figure [l] where we show the strain 
amplitude ratio of higher modes relative to the (2,2) 
mode. The top panel shows the case of an equal mass 



BBH and the bottom one with a 1:4 mass ratio. It is 
evident in both cases the dominance of the (2,2) mode. 
However, for the unequal mass case, the second strongest 
mode, the (3,3) mode, is more than 10% of the (2,2), 
reaching above 20% near merger. On the other hand, 
for the equal mass case, the second strongest mode is 
the (4,4) mode, and it is only a few percent of the (2,2) 
in strength. Similar situations occur for precessing (i.e. 
spinning BHs) binaries. 

If source orientation effects are considered, the relative 
strength of modes can also be significantly modified. Re- 
cently, we carried out a study that addresses such effects 
for quadrupolar templates [2J. We found that for pre- 
cessing binaries, overlaps could drop below 0.97 for up to 
65% of the sky at the source location. The present study 
takes the next step. We investigate overcoming the lim- 
itations of quadrupolar templates by identify hierarchies 
of higher modes that could be added to the templates, 
increasing their sky coverage. Our study assumes that 
the detector is optimally oriented, and we focus on line- 
of-sight, source orientation effects. Therefore, "the sky" 
coverages reported will be from a source-centric frame 
of reference. Furthermore, we focus on BBHs with total 
mass M = 1OOM0. For these relatively massive binary 
systems, interferometers such as LIGO only "see" the late 
inspiral and merger of the binary. Therefore, templates 
in this study are obtained entirely from NR simulations. 
For less massive binaries, the early inspiral becomes more 
important. The templates in those cases would have to 
be obtained from hybrid waveforms that stitch NR and 
post-Newtonian strains [3l H] . 

Our study finds that template banks for GW detec- 
tion based only on the (2,2) mode are probably suffi- 
cient for comparable mass, low precession BBHs. At the 
same time, our work demonstrates that the effectiveness 
of match filtering is severely impaired for highly precess- 
ing or moderately large unequal mass BH systems if tem- 
plates banks are based entirely on the quadrupolar mode. 
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FIG. 1: Strain mode amplitude ratios relative to the (2,2) 
mode non-spinning BBHs. Top panel shows the case of an 
equal mass binary, and the bottom panel that of a 1:4 mass 
ratio system. Modes not included have relative amplitudes 
less than 10" 3 . 



II. MATCHED FILTERING AND DETECTION 
REACH 

A GW signal s impinging an interferometer has the 
following structure 

S (t;6, 0,^,1) = F+(9j^)h + (t;0 

+ F x (0,$,rl>)h x {t;i), (1) 

where h + and h x denote the two polarization strains of 
general relativity. h + and h x depend on a parameter 
vector £. Some of the components in £ are parameters 
intrinsic to the system, such as the binary eccentricity, 
BH masses and spins; others are external parameters like 
the binary orientation and its distance to the detector. 
In addition, the signal s depends on the characteristics 
of the detector through the response or antenna pattern 
functions F + and F x jS] . For an interferometer with arms 



oriented along the x and y axes, these functions read 

F + = — -(1 + cos 2 9) cos 2<f> cos 2if> 

— cos 9 sin 2(f> sin 2ip (2) 
F x = ^(1 + cos 2 0) cos 20 sin 2^ 

— cos 9 sin 2(f> cos 2ip , (3) 

where the angles (8, </>) are the location of the binary sys- 
tem in the sky of the detector and %f> the polarization 
angle. A detector is optimally oriented if a signal arrives 
with 9 = or it and = 0. On the other hand, the 
detector is blind for signals with 9 = tt/2, arriving from 
directions aligned with the arms of the interferometer. 
With a network of interferometers, the issue of detector 
blindness can be alleviated [SJ [7] . For this reason, our fo- 
cus will be on "intrinsic" source orientation effects, which 
depend only on the direction of propagation of the GW. 
That is, the signals arriving at the detector will be 

s(t;8,<f>,r,£ b ) = h + (t;9,(/),r,£ b ) cos2<0 

+ h x (t;8,(f>,r,Cb) sin2V>. (4) 

In Eq. Q, the angles (9, </>) give the direction of propaga- 
tion of the GW pointing to the detector from a reference 
system located at the source, r is the distance between 
the source and the detector, and the vector denote the 
parameters intrinsic to the binary (spins, masses and ec- 
centricity). For the present study, the reference system 
used at the source is the one used to carry out the NR 
simulations. That is, the reference system has its origin 
at the center of mass of the binary, and has its z-axis 
aligned with the orbital angular momentum at the start 
of the simulation. 

Given the Fourier transforms A and B of two real time- 
dependent signals or waveforms, the inner product with 
respect to the noise spectrum density S n (f) of the detec- 
tor is given by 

With this inner product, the signal-to-noise (SNR) of a 
signal s and a template u is given by [5] 

, . ( s \ u ) /„s 

p{s, u) = max — === . (6) 

tO.V'O y/(u\u) 

The maximization over to enters because the signal ar- 
rives at an unknown time. Thus, one needs to "slide" the 
template relative to the signal, which in Fourier space im- 
plies replacing u(f) by u(f) exp(— 27ri/i ) in Eq. In 
principle, there is also an unknown phase difference ipo 
between the signal and the template related to the po- 
larization angle ip. This introduces an additional factor 
of exp(2iriif»o) . If u = s in one gets the so-called 

optimal SNR p(s, s) = yj (s|s). 
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The match or overlap p between the signal s and a 
template u is given by 



p(s,u) 



(s\u) 



to, A> y/ (s\s) (u\u) 



(7) 



Thus, p(s, u) — p(s, w) y/ (s\s) — p(s, u) p(s, s). In some 
instances, it is more convenient to work with mismatches, 
which are defined as e = Notice that the 1/r scaling 

of the signal s implies that p oc 1/r. However, because 
of the normalization factors, p and e are independent of 
r. Furthermore, the distance scale r used to construct a 
template does not play any role because of the normal- 
ization factor yj (u\u) in Eqs. ^ and ^7j. 

Consider two signals s\ and S2 from two identical GW 
sources, S\ located at distances r\ and S2 at r 2 - Then, 



p(si,u) _ p{si,u) _ pjsi^Sj) 

p(s 2l u) p(s 2 ,u) p(s 2 ,s 2 ) /o(s 2 ,s 2 )' 



(8) 



where the last equality follows because the sources are 
identical, and thus p(si,u) — p(s 2 ,u). On the other 
hand, 

p(s 2 ,s 2 ) = —p(r 2 s 2 ,r 2 s 2 ) 
r 2 

l n 

= —p(r 1 s 1 ,r 1 s 1 ) = —p(s 1 ,s 1 ), (9) 
r 2 r 2 

where we have used that the optimal SNR of a signal 
times its distance, r s, is independent of r. Therefore, we 
can rewrite Eq. (pi) as 



p(si,u) _ p{sy 1 s 1 ) _ r2 
p(s 2 ,u) p(s 2 ,s 2 ) n 



(10) 



It is wrong to conclude from Eq. ( 10 1 that the ratio r 2 /ri 



depends on the template u because of p(s\,u)/ p(s 2 ,u). 
The dependence on u is eliminated because we are consid- 
ering the same source, just located at different distances. 

The SNR property fcty allow one to define for a signal 
s a horizon distance R[s) as follows: 



R(s) 



p(s, s) 
p(s,s) 



(11) 



where p(s, s) is the optimal SNR of the "signal" s = r s. 
R(s) should be interpreted as providing the maximum 
distance that an interferometer is able to detect a signal 
s given an optimal SNR threshold. For instance, during 
LIGO's science run S6, the horizon distance of a com- 
pact binary coalescence with optimal SNR threshold of 
p(s, s) = 8 Mpc -1 was estimated to be R(s) s» 40 Mpc 
using as a model for the signal a post-Newtonian strain in 
the stationary phase approximation [9]. Therefore from 
Eq. Jul), p(s,s) « 320. 



For the present study, we are interested in investigating 
the effect that the choice of a template u has on the abil- 
ity of an interferometer to detect a signal s. Therefore, 



instead of the horizon distance R(s), which is indepen- 
dent of the template, we need a template reach distance 
lZ(s,u). Since the mismatch p characterizes the "prox- 
imity" of a template to a signal, we define the template 
reach distance as 



1Z(s, u) = p(s, u) R(s) . 



(12) 



Given that p(s,u) — p(s,u), we can rewrite Eq. (12) in 
a form similar to that of Eq. ( [TT| ) : 



1Z(s, u) 



p{s,u) 
p(s,s) ' 



(13) 



Notice that TZ(s, s) — R(s) if u = s; otherwise, TZ(s, u) < 
R(s). Furthermore, the relative change in reach is given 
by the mismatch: 



6R _ R(s)-K(s,u) 



R 



R(s) 



1 — p(s, u) = e(s, u) . (14) 



III. BINARY BLACK HOLE WAVEFORMS 

As mentioned in the introduction, the BBHs we con- 
sider have total mass M — 100 Mq. For these binaries, 
interferometers such as advanced LIGO are most sen- 
sitive to GW frequencies emitted by the binary during 
the late inspiral and merger; and, therefore, the tem- 
plates and signals needed in our study are those con- 
structed entirely from NR simulations. The simulations 
were produced with the Maya code [TUHT5] of the NR 
group at Georgia Tech. Maya uses the Einstein Toolkit 
[To] which is based on the CACTUS [T7] infrastructure and 
CARPET [TS] mesh refinement. Evolution thorns were gen- 
erated with the Kranc code generator. 

We considered a variety of BBH configurations in 
quasi-circular orbits. Table [TJ gives the characteristics of 
each binary: BH mass ratio q = mi/m 2 , dimensionless 
spin parameters x, angle C between the BH spin and the 
z-axis in the xz-plane, initial separations d/M, and grid 
spacing A/M in the finest refinement level, and the cov- 
ering factors r 2 .2 when only the (2,2) mode is used as 
the template and r 2i * which includes all the \m\ < I — 2 
modes in the template (see next section) . The configura- 
tions are classified into three series. Series Q consists of 
BBH with non-spinning BHs. For the series S, at least 
one of the BHs is spinning and aligned with the orbital 
angular momentum; thus, there is no precession. Finally, 
the series P consists of precessing binary systems. 

The templates and signals were built from strains /i+ 
and h x obtained from the Weyl Scalar ^4, one of the 
outputs of the simulations. Specifically, we decompose 
the Weyl Scalar \l/4 produced during the simulation into 
spin- weighted spherical harmonics as 



l.m 



(15) 
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0.6 
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0.247 
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0.6 
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10 
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4.0 


0.6 
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0.188 


0.234 


P22 


4.0 


0.6 





6 
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10 
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0.130 


0.218 


P23 


4.0 


0.6 





6 


270 


270 


10 


140 


0.026 


0.244 



TABLE I: Characteristics of the BBH systems used in this 
study, q = m 1 /m2 denotes the BH mass ratio, Xi,2 the di- 
mensionless spin parameter, £1,2 the angle between the BH 
spin and the z-axis in the zz-plane, d/M the initial binary 
separation, and A/M the grid spacing at the finest mesh. 
The last two columns give the covering factor T2,2 when only 
the (2,2) mode is used as the template and the covering factor 
T2,* when the template includes all the jmj < I — 2 modes. 



where the angles 6 and <p are relative to a coordinate 
system with origin at the center-of-mass of the binary and 
with a z-axis aligned with its orbital angular momentum 
at the beginning of the simulation. Then, we solve ^4 = 
h + — i h x = h* to get the strains 



M 



h(t;6,cf>) 



E 



(16) 



For our signals s, we calculate h + containing all the 
(£, m) modes with 2 < £ < 5. The templates u are con- 
structed from different subsets among all the 32 (£, m) 
modes for the signals. Explicitly, a template u n contain- 
ing n < 32 modes is given by 



4 = 1 



(17) 



with 
i.e., 

i collectively labels (£, m) 



Wi denoting a mode in the sense of Eqs. (16 1, 
Wi = Hi m {t) - 2 Y l l k m (9 7 (j))M/r in which the index 



To avoid the integration to 
strain, and the errors associated with it, we calculate the 
strain in the frequency domain for the matches, that is, 
h = ^/(-4tt 2 / 2 ) • 



IV. COVERING FACTORS AND MODE 
HIERARCHIES 

As mentioned in Sec. [nj the focus of our analysis is the 
mismatch e of a signal s(£ s ) with a template m(£ u ), where 
the parameter vectors £ s and £ u include both, intrinsic 
physical (BH masses and spins) and extrinsic (position 
of the source in the sky of the detector and orientation 
of the binary) parameters. Since we are only concerned 
with line-of-sight orientation effects, we will assume that 
the detector is optimally oriented and £ s = £ M , and in 
particular (8 s ,(j) s ) = (8 

ui4*u)- As a consequence, once 
a source is selected, e will only depend on the location 
of the source in the sky, i.e. e = e((9, <f) as implied by 
Eq. @. 

We will mostly present results as sky-maps, and inves- 
tigate how e(9, <fi) depends on the particular content of 
higher modes in the template. We construct sky-maps 
with 61 x 61 pixels; each pixel "colored" with its value 
of e. Figure [2] gives an example of two sky-maps of mis- 
matches for a BBH with mass ratio q = 4. The top panel 
shows a non-precessing binary (case Q06 in Table [TJ) m 
which the template contains only the (2,2) mode. The 
bottom panel shows a precessing binary (case P19 in Ta- 
ble m. The template in this case contains all the m- modes 
for £ = 2. Notice that, the mismatch is smaller around 
the poles in the non-precessing case. On the other hand, 
for the precessing case, the region where the mismatches 
are smaller have been shifted due to the precession of the 
system. 

In addition to the sky-maps, given a mismatch toler- 
ance e mx , we label a pixel on if e < e mx and off other- 
wise. With the on-pixels, we estimate the covering factor 
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FIG. 2: Representative sky-maps of mismatches for a BBH 
with mass ratio q = mi /m,2 = 4. The top panel shows a non- 
precessing binary (case Q06) in which the template contains 
only the (2,2) mode. In the bottom panel, the precessing case 
P19 is shown using a template that contains all the m-modes 
for i = 2. 



the value t 2 . This process is repeated until we obtain a 
template, un, such that tn > T tg . 

Second Stage: Because the procedure to get {u n } is 
locally optimal, it is not guaranteed that the sequence 
has the minimum possible number of steps to reach Ttg, 
in other words, that the template un has the minimum 
number of modes. Next is to identify all modes in un 
that could be excluded and still have tn > T tg . Starting 
with un, we construct templates 



UN 



N 



N 



(18) 



for n = 1 . . . N. We then select the it* for which t* > 
T tg . If there is more than one template it* , we single 
out the case in which the highest mode was subtracted. 
That mode is then remove from the template un and 
set ujv-i = u* n , which contains N — 1 modes. The same 
process is applied now to un-i- That is, we construct 
templates 



T^mx) of a template as the fraction of the sky where 
£ < £mx- This means that a template that cannot de- 
tect any signals would have a covering factor of 0, and 
a template that matches all signals would have a cov- 
ering factor of 1. From now on, when stating covering 
factors, we will just use t, understanding that e mx is im- 
plied. Unless explicitly specified, we use e mx — 0.03 in 
the remainder of the paper, a customary value within the 
GW data analysis community. For instance, the covering 
factors for the cases in Figure [2] are t = 0.177 (top) and 
t = 0.242 (bottom). 

It should not be surprising that the mode content in a 
template has an impact on its mismatch with the signal, 
and thus its covering factor. Our goal is then to find the 
optimal template un that yields a target covering factor 
Ttg. We use a greedy algorithm for this purpose. The end 
result is a sequence of templates {u n } built from a set of 
modes {w n } as in Eq. (17), with n = 1... N. At each 
step in the sequence, a template u n is locally optimal 
in its ability for improving sky coverage and reaching 
Ttg. Furthermore, the template sequence is ordered in a 
hierarchy with respect to their sky coverage capabilities. 
The procedure we use to obtain this hierarchy involves 
three stages. 

First Stage: The first template u\ is selected to be the 
dominant (2,2) mode. The covering factor for this first 
template is denoted by T\. The second template u 2 is 
obtained from a superposition of the previous template, 
ui, and a mode Wi selected from the remaining 31 modes. 
This yields 31 potential candidates for u%. We promote 
to template it 2 the candidate that delivers the best im- 
provement in sky coverage over the value t\. The third 
template in the hierarchy is constructed in a similar fash- 
ion. That is, 1*3 = u 2 plus one of the 30 remaining Wi 
modes. Among the 30 candidates for 113, we keep the 
one with the largest improvement in sky coverage over 



UN-l - w n 



N-1 

£ 

i=l 



N-1 

£" 



(19) 



for n = 1 . . . N— 1. The outcome is a new template v-n-2- 
The process is repeated L times until we are not able to 
construct templates it* such that t* > T tg . The resulting 
u m j with M = N — L, is the template with the minimum 
number of modes having a covering factor T tg . 

Third Stage: With the template Um a t hand, we ap- 
ply a greedy algorithm to re-construct the hierarchy but 
now in reverse order; that is, starting with um we obtain 
um-i as the template that decreases f M the least. We 
continue until we reach a template made out of a single 
mode. 

Tables [TTJ |III| and IV show how our greedy algorithm 
builds a hierarchy for the P17, P19 and P20 cases, re- 
spectively. The top row labels the (£, m) modes, and the 
left column denotes the covering factor t. The target 
covering factor is r tg — 1. For a given r, the Xs in the 
row indicated the modes that were included in the tem- 
plate. The Os denote the mode from the previous step 
that was removed. The first horizontal line denotes the 
end of the first stage in the greedy algorithm. The rows 
between the first and second line are the steps involved 
in the second stage. Finally, below the second line are 
the templates in hierarchy one is seeking, in decreasing 
order of sky coverage. 



V. RESULTS 

A sense of the effectiveness of the mode (2,2) in cov- 
ering the sky is given in the second to last column of 
Table [I] where we report the covering factor of just the 
(2,2) mode, given by t 2j 2- As long as q < 1.5, the (2,2) 
mode will effectively cover the entire sky for non-spinning 
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TABLE II: P17 Hierarchy: Table showing how our method for constructing the mode hierarchy works for the case P17. The 
top row labels the (£, m) modes, and the left column denotes the covering factor r. The target covering factor is T tg = 1. For a 
given t, the Xs in the row indicated the modes that were included in the template. The Os denote the mode from the previous 
step that was removed. The first horizontal line denotes the end of the first stage in the greedy algorithm. The rows between 
the first and second line are the steps involved in the second stage. Finally, below the second line are the templates in hierarchy 
one is seeking, in decreasing order of sky coverage. 
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TABLE III: P19 Hierarchy: Same as in Table [TT] but for the P19 case. 



and non-precessing BHs (see Q and S series). Even in the 
case of moderate precession (e.g. cases POland P02) the 
quadruple mode is sufficient. At the same time, it is evi- 
dent from the Q-series that as q increases the effectiveness 
of a (2,2)-only template quickly deteriorates. For q = 4, 
the covering factor is already r = 0.177, reaching 0.097 
for q — 10. In the last column of Table |TJ we report 
the covering factor T2 ; * when the template includes all 



|m| < I = 2 modes. Notice that for the non-precessing 
binaries (series Q and S), the inclusion of these modes 
does not translate into a significant improvement of sky 
coverage. For the precessing systems, ti^i is not a consis- 
tent description of the coverage, because it depends on 
the fixed extraction frame. A frame independent choice 
is T2,*. Typically, T2,* improves coverage by at least a 
factor of 2 for the precessing cases, and in the highly 
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TABLE IV: P20 Hierarchy: Same as in Table [TT] but for the P20 case. 



precessing cases, by a significant amount more. 



A. Q-series 

We first analyze the impact of adding other modes to 
the quadrupolar templates for the Q-series. To anticipate 
what to expect from the greedy algorithm described in 



Sec. IV we will consider first adding modes based only 
on their energy radiated; these are (3,3), (4,4) and (2,1) 
for the Q-series. Adding the (3,3) mode translates into 
a significant improvement as depicted in the top panel 
in Figure [3] where we show a bar chart of the covering 
factor as a function of the mass ratio. For each bar, black 
denotes the covering factor that is achieved with only the 

(2.2) mode. The gray bar above black is the enhancement 
that one gets by considering a template with (2,2) + 

(3.3) . One sees that now full coverage is obtained for 
q < 3 and substantial coverage improvement of 0.5 < 
t < 0.75 for 7 > q > 5, and r = 0.45 for q = 10. Adding 



the next most energetic mode, the (4,4) mode, produces 
coverage above 0.8 for the remaining cases except q = 10, 
as seen by the white bars in Figure [3] top panel. Full 
coverage for all the mass ratios considered is achieved 
when the (2,1) is added (light gray bars). 

The ordering based on energetic used for the Q-series 
gets slightly modified if one uses our greedy algorithm. 
The method reveals that adding the (2,1) mode before 
the (4,4) mode is more optimal in reaching full cover- 
age, as seen in the bottom panel of Figure [3j The rea- 
son why the greedy algorithm favors the (2,1) mode for 
q > 4 is because, although the mode is less energetic, it 
channels energy into regions of the sky not covered by 
the (2,2)+(3,3) template more effectively than the (4,4) 
mode. This can be seen in Figure [4] for the q — 10 case in 
the Q-series. Each of the panels in this figure shows the 
percentage of total energy emitted in a given direction 
that is channeled through the modes present in a tem- 
plate. From top to bottom, the panels show the cases 
for templates (2,2), (2,2)+(3,3), (2,2)+(3,3) + (4,4), and 
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FIG. 3: Covering factors as a function of mass ratio for the 
Q-series. Each portion of the bar chart represents the con- 
tribution to the covering factor by adding a mode. In each 
bar, bottom to top is the order used in adding modes. The 
gray scale is as follows: black (2,2) mode, gray (3,3) mode, 
white (4,4) mode and light gray (2,1) mode. Top panel shows 
the case in which modes are added according to their energy 
emission strength. The bottom panel shows the hierarchy or- 
der obtained by our greedy algorithm, which translates into 
reversing the order of the (4,4) and (2,1) modes. 



(2,2) + (3,3)+(2,l), respectively. The black lines in each 
figure denote mismatches e = 0.03 between the template 
and the signal. Darker regions have mismatches e > 0.03. 
Thus, not surprisingly, there is a correlation between the 
amount of energy captured by the template and its mis- 
match value. Roughly speaking, regions for which the 
template captures more than 80% of the energy emit- 
ted in a given direction have e < 0.03. In these regions 
of low mismatch, the average percentage radiated in all 
four panels of Figure|4]is approximately 90% overall. It is 
then evident from the bottom two panels in Figure|4j that 
adding the (2,1) mode before the (4,4) yields a larger re- 
gion where the template better captures the energy emit- 
ted. 

We have also found that there is a degree of degen- 
eracy. For instance, in the q — 7 case, adding to the 




FIG. 4: Percentage of total energy emitted in a given direc- 
tion that is channeled through the modes present in a tem- 
plate for the q = 10 case in the Q-series. From top to bottom, 
the panels show the cases for templates (2,2), (2,2) + (3,3), 
(2,2) + (3,3) + (4,4), and (2,2) + (3,3) + (2,l), respectively. The 
black lines in each figure denote mismatches e = 0.03 between 
the template and the signal. Light regions with values larger 
than 0.8 have mismatches e < 0.03. 

template (2,2)+(3,3)+(2,l) the mode (4,4) or the mode 
(5,5) accomplishes full coverage. That is, both modes 
contain enough power in comparable parts of the sky. 



B. S-series 

Next to analyze is the S-series, consisting of BBHs with 
spinning BHs aligned with the orbital angular momen- 
tum, i.e. non-precessing binaries. As mentioned before, 
for low mass ratio q < 1.5, the quadrupolar mode is able 
to get good coverage of the sky. We will focus then on 
the q = 4, 5 and 6 cases in Table [TJ These are binaries 
with only one of the BHs spinning. As with Figure [3j the 
bar charts in Figure [5] show the incremental effectiveness 
of adding modes to a template. In this case, the cover- 
ing factor is plotted as a function of the dimensionless 
spin parameter \+- The gray scale is: black (2,2) mode, 
gray (3,3) mode and light gray (2,1) mode. From top 
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to bottom, the panels show the q = 4, 5 and 6 cases, 
respectively. 

The first thing to notice is that if only the (2,2) mode 
is used as a template, the covering factors (black bars) 
are basically independent of the spin. In addition, the 
values of the covering factors for the (2,2) mode template 
decrease as the mass ratio increases. The big changes 
occur when one adds the (3,3) mode to the quadrupolar 
template. In some instances, one is able to increase the 
coverage to include the entire sky. It is also interesting 
that the influence of this mode increases with the value of 
the spin, consistent with the energy of the (2,2) and (3,3) 
modes increasing with spin. For example, in the 9 = 4 
case, a (2,2)+(3,3) template is able to saturate the sky 
for x+ > 0-4- Finally, for those cases in which the (3,3) 
is not able to yield full sky coverage, i.e. low spin, high 
mass ratio binaries, adding the (2,1) mode completes the 
sky. 



The last series focused on precessing binaries. In this 
series, our greedy algorithm used to construct hierarchies 
provided insights on the number of modes that will be 
needed for detecting GW from astrophysically realistic 
binaries. We applied the greedy algorithm to each of the 
23 cases in the series. However, we show the steps that 
the algorithm takes for only the P17, P19 and P20 cases. 
The steps are shown in Tables |Tl} III and IV respectively. 



A dramatic aspect found for these highly precessing 
binaries is that it takes a substantial number of modes 
to cover at least 60% of the sky. For the P17 case, 
one sees from Table [TT] that one needs a template with 
8 modes: (2, ±2), (2,±1), (2,0), (3,3), (3,-2) and (3,0). 
Similarly, for the P19 in Table [111] the template contains 
the (2, ±2), (2,±1), (2,0), (3,3) and (3,0) modes. Finally, 
for the P20 case, as seen from Table |HI[ the template 
involves the (2, ±2) and (3, ±3) modes. The entire hier- 
archy for each case is listed below the second line in each 
of the Tables ordered in decreasing covering factor. It 
is clear then that precession plays a big role on the out- 
come of the hierarchy, and in particular the modes that 
are needed for full coverage of the sky (i.e. the row below 
the second line in the tables). 

To quantify the influence of precession, Figure [6] shows 
the total number of modes needed for full coverage in the 
P-series organized by mass ratio value q and spin \. The 
number of modes are given as a function of the angle 
where cos© = (L ■ S)/(LS). Here L is the orbital 
angular momentum and S — S\ + S2 the total spin at 
the beginning of the simulations. Between 0° < < 90°, 
the larger the angle the higher the observed precession 
in the system; and, therefore, it is not surprising in this 
range of angles for Figure[6]to indicate that the number of 
modes needed for full sky coverage grows monotonically 
with 0. One can also see that, for a given angle value, 
the number of modes needed also increases monotonically 






FIG. 5: Covering factors versus spin, \+ = a +/ m +i of the 
more massive BH. The less massive BH has vanishing spin. 
The top panel is q — 4, center q = 5, and bottom q = 6. The 
gray scale is as follows: black (2,2) mode, gray (3,3) mode, 
and light gray (2,1) mode. 



with the mass ratio q. For angles > 180°, we only have 
a few data points. But it is clear that the monotonicity no 
longer applies since for those cases the level of precession 
decreases. 

To check the robustness of our algorithm in finding 
hierarchies, we selected two q — 4 binaries, one non- 
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FIG. 6: Optimal number of modes needed for full coverage 
versus O = acos^L ■ S/(L S)) for all processing simulations 
in the P-series. 

precessing (Q06) and the other with one of the high- 
est precession (P23). We ran these two cases at dif- 
ferent resolutions: (M/160, M/180, M/200) for Q06 and 
(A//140, M/160, M/180) for P23. In addition, we con- 
structed the GWs from the simulations at three different 
extraction radii (r = 60, 75, 90M). Our greedy algorithm 
was able to follow identical steps and find the same hi- 
erarchy for the Q06 case. For the precessing case P23, 
there were small differences in the steps taken for the 
lowest resolution and smallest extraction radii case, but 
the resulting hierarchy remained the same. 

We also experimented with steps that add and remove 
(£, \m\) and (£, — \m\) modes simultaneously. This re- 
duced the number of total possible steps from 31 to 17 
after the initial seed of (2, ±2). The net effect was to 
shrink or eliminate Stage 2. We also found that there 
were very few cases where the resulting hierarchy con- 
tained fewer modes than the single mode case. 



VI. CONCLUSIONS 

As the NR community explores the complete, generic 
parameter space of BBH spacetimes, non-dominant 
modes play an increasingly important role in detect- 
ing and characterizing potential gravitational wave sig- 
nals. In this paper, we introduced a hierarchical, greedy 
method to identify the modes necessary for full sky cov- 
erage. We define sky coverage as the percent of sky in 
a source-centric frame that a template will match with 
nature up to a mismatch threshold, here taken as 0.03. 

The loss of sky coverage is most noticeable in highly 
unequal and precessing BBH systems. As long as q < 
1.5, the (2,2) mode will effectively cover the entire sky 
for spinning and mildly precessing systems. For non- 
spinning systems, as the mass ratio increases to q > 7, 
the inclusion of four modes is necessary. In the absence 



of precession, a larger spin magnitude increases the en- 
ergy of the system and fewer modes are necessary. As 
precession becomes more dramatic, as measured by 0, 
the range of sky covered decreases to below 60-90% and 
the number of modes necessary to achieve full coverage 
increases to 12 for the many of the cases studied here. 
For non-precessing systems, the four important modes 
are also the four most energetic: (2,2), (3,3), (4,4), and 
(2,1). For all precessing systems covered here, all the 
£ = 2 and £ = 3 modes are sufficient for covering all 
possible orientations. 

Although we have focused on binaries with total mass 
of 100 M , our methodology is directly applicable to less 
massive systems, where templates are constructed from 
hybrid waveforms that combine post-Newtonian and NR 
results into one, long waveform. In broad terms, our con- 
clusions regarding the influence of higher modes are likely 
to carry over to lighter binaries, see for instance Fig- 
ure 1 and related work in [20 j . and thus have an impor- 
tant impact for the enterprise of building template mod- 
els. Our method could potentially give clues where effort 
in analytic modeling, NR simulations, and hybridization 
techniques will yield the most benefit. In particular, it 
should be possible to estimate the loss in sky coverage 
due to missing information resulting from using a subset 
of modes. 

In current GW data analysis searches based on 
matched filtering, waveform models are used to construct 
a bank of templates placed at discrete points in the two- 
dimensional parameter space spanned by total mass M 
and mass ratio q |21j . Such a bank responds not only 
to GW signals, but also to noise in the detector, cre- 
ating a population of false alarm triggers in the search 
pipelines. As a consequence, a true signal must stand out 
sufficiently above those triggers in order to be claimed 
as a detection. If one were to utilize higher-mode in- 
formation, it would become necessary to increase the 
dimensionality of the template bank by two to include 
the orientation parameters 9 and <fi. There is a poten- 
tial tension in doing so. Although the additional tem- 
plates will assist in responding better to the presence of 
signals, the increased number of templates will unfortu- 
nately also elevate the population of false alarm events 
thus impairing detection ability. In order to maximize 
the template response and minimize the emergence of 
false alarm triggers, one would have to adjust various as- 
pects in the current search methodologies. For example, 
currently two distinct searches are performed, one cov- 
ering the mass region M < 25Af Q and another spanning 
25M < M < IOOMq. This is done in part because the 
background populations generated by templates in these 
two regions differ significantly. It is possible that higher- 
mode information could best be employed by further sub- 
dividing by mass ratio. The NINJA-2 project [22 , which 
is currently ongoing, provides one possible route to in- 
vestigating such issues. NINJA-2 will be able to add 
complete hybrid waveforms which include higher modes 
into realistic next-generation detector noise. Using the 
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resulting combined data with guidance from techniques XSEDE TG-PHY120016 and the Cygnus cluster at Geor- 
such as the one introduced here may be able to help as- gia Tech. 
sess the difficulties and rewards of higher- mode templates 
in a complete search. 
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